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Abstract 

We report about recent progress in the treatment of bound states in QCD. 
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on 1 Introduction 

One of the most difficult and less understood problems in Quantum Field Theory is 
the treatment of bound state systems. Even in QED where the interaction is simply 
given in perturbation theory, the practical calculation of the bound state properties is 
tricky The complication comes from the mixing of the many characteristic scales of 
the bound state, the mass of the particles, the momentum and the energy of the state. 
A way to handle this problem is provided by Non Relativistic QED (NRQED) which 
is an effective theory equivalent to QED and obtained from QED by integrating out 



the hard energy scale m. The ultraviolet regime of QED (at energy scale m) is 



ON 
O 



perturbatively encoded order by order in the coupling constant a in the matching 



coefficients that appear in front of the operators in the effective Lagrangian. Each 
term in the effective Lagrangian has a definite power counting in a and then a 
disentangling of the various scales is achieved. A similar simplification is obtained 
using Non Relativistic QCD (NRQCD) in heavy quark systems. These systems 
are characterized by a dynamical dimensionless parameter, the quark velocity v, 
^ . which is small and allows a classification of the energy scales of the problem in 

hard (~ m), soft (~ mv) and ultrasoft (~ mv 2 ). This provides a power counting 
scheme. Here, however, a further and conceptual difficulty arises in connection 
with the nonperturbative nature of low-energy QCD. The relation between v and 
the QCD parameters is unknown but certainly v includes both perturbative and 
nonperturbative effects. Therefore, the evaluation of the bound state properties is 
ultimately done with a lattice simulation. In the case of bound systems with at 
least one light quark simplifications of this type do not hold. On one hand no small 
expansion parameter exists, on the other hand the light quark mass is generated 
via chiral symmetry breaking. Hence, the interplay between confinement and chiral 
symmetry breaking has to be considered. 

In this talk we report about recent progress in the treatment of bound states in 
QCD. We show that it is possible to obtain a model-independent and gauge-invariant 
result for the heavy quark interaction at order t> 4 of the systematic expansion in v. 
The interaction turns out to be simply given in terms of a generalized (distorted) 
Wilson loop. The result is suitable for lattice evaluation as well as for analytic 
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evaluation once a QCD vacuum model is considered. We show that the results for 
the heavy quark dynamics are substantially under control and are given in terms of 
two nonperturbative parameters T g , the gluon correlation length, and G 2 , the gluon 
condensate. Adopting the same framework in order to study the heavy-light bound 
states in the non-recoil limit, spontaneous chiral symmetry breaking and a confining 
chiral non-invariant interaction emerge quite naturally. We discuss this last case 
with more details. 



2 The heavy quark interaction 

The heavy quark interaction can be obtained analytically at the order v 4 of the 
systematic expansion of the interaction in v. Here, we report only the main steps of 
the derivation referring to |]1], |2| for further details. 

• Step 1. Set up the NRQCD Lagrangian. The NRQCD Lagrangian [[| is 
obtained from the QCD Lagrangian via a Foldy-Wouthysen transformation. 
At order 0{v A ) the NRQCD Lagrangian describing a bound state between a 
quark of mass mi and an antiquark of mass m 2 is [|J] 



2mi 8mf * 2m 1 u 8mf 
+zcl 1) ^ (T ^ X "^ 9 ^ X ^ }Qi + antiquark terms (1 <-> 2) H — — Q1Q2Q2Q1 

8mf J 171x1712 

+-^—Q\aQ 2 QlaQ 1 + -^—Q\T a Q 2 QlT a Q 1 + -^Q\T a aQ 2 Q\T a aQ 1 , 
mim 2 m\m,2 mii7i2 

where Qj are the heavy quark fields. Reparameterization invariance[|J fixes 
c 2 = c 4 = 1. The coefficients Cp\ , ... are evaluated at the matching scale fi 
for a particle of mass mj. They encode the ultraviolet regime of QCD order by 
order in a s . The explicit expressions and a numerical discussion can be found 
in 0. The power counting rules for the operators of Eq. (jH) are Q ~ (mv ) 3//2 , 
D ~ mv, gAo ~ mv 2 , gA ~ mi; 3 , ~ m 2 v 3 and (/I? ~ m 2 v A . Four 
quark operators which are apparently of order v 3 are actually suppressed by 
additional powers in a s in the matching coefficients and the octet contributions 
by an additional power in v 2 on singlet states. Therefore in the following we 
will neglect these contributions with the exception of a term which mixes 
under RG transformation with the chromomagnetic operator contribution to 
the spin-spin potential 0. We will call the corresponding matching coefficient 
d. 

Step 2. Write down the quark- antiquark gauge-invariant Green function using 
a path integral representation. The 4-point gauge invariant Green function G 
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Figure 1: Distorted Wilson loop. 

associated with the Lagrangian (P is defined as [7]. |[ 

G(x 1 ,y 1 ,x 2 ,y 2 ) = (0\Ql(x 2 )U(x2,x 1 )Q 1 (x 1 )Ql(y 1 )U(y 1 ,y 2 )Q2(y2)\0), 

where U(x 2 ,xi) = exp l—ig J ds (x 2 — x\YA^{xi + s(x 2 — ^i)) j is a Schwi- 

nger line added to ensure gauge invariance. After integrating out the heavy 
quark fields, G can be expressed as a quantum-mechanical path integral over 
the quark trajectories: 



Zj - rrij - ' 
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G(x 1 ,y u x 2l y 2 ) = / Vz{Dp x I Vz 2 Vp 2 exp < i / dtJ]Pj 

J Vi J V2 [ J -T/2 j=1 
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x exp [ T/2 dtg 2 dT*»*WT a ®*®5 3 (z 1 - z 2 

I m\m 2 J-T/2 
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where the bracket means the Yang-Mills average over the gauge fields, T is the 
Wilson loop made up by the quark trajectories Z\ and z 2 and the endpoints 
Schwinger strings and y\ — y\ = —T/2, x^ = x\ = T/2, see Fig. [I| 

S'tep 3. Extract the form of the heavy quark interaction from Eq. (0). Assum- 
ing that the limit exists, we define the heavy quark- ant iquark potential V as 

rT/2 

lim / dtU/T. Expanding in v we get 

T->co J-T/2 

V(r) = lim "° g(ly(r)) + (3) 

v ' r-00 T 

/SWJLW S( 2 ) -L( 2 ) \ 24^/(^ + 4^(0 ■ L( 2 ) + ■ LW 4^(r) 

V m 2 m 2 / 2r m\m 2 r 
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where iy(r) = P exp |— ig j> dz^A^z)^ is the averaged value of the de- 
formed Wilson loop, see Fig. []]. The expansion of it around the static Wil- 
son loop ^(To) (r is a r x T rectangle) gives the static potential Vo = 
lim ilog(W(r Q ))/T plus velocity (non-spin) dependent terms [ED, |9| which are 

controlled by four scale- independent potentials Vj(r) i = b,c, d, e. and 
are the spin and orbital angular momentum operators of the particle j. The 
matching coefficients are defined as 2dp S = Cp S ± Cp S . The spin-dependent 
potentials agree with the ones obtained in refs. || with the exception of the 
matching coefficients that were introduced in [H]. 

All the potentials Vi — V4 V a — are obtained as functions of the average value of 
the distorted Wilson loop (W(T)) and insertions of one or two field strengths in this 
average, (F^WiT)) or {F^Fp^W (Y)) . Via deformation of the quark or (antiquark) 
trajectory, such v.e.v. of field strength insertion can be calculated via functional 
derivatives of the Wilson loop. We conclude that to obtain the complete quark- 
antiquark order 0(v A ) interaction (quenched) no other assumptions are needed than 
the behavior of (W(T)): given (W(T)) everything is analytically calculable. On 
the other side, expanding the average of the distorted Wilson loop T in terms of 
the static Wilson loop T we get expressions for the potentials suitable for lattice 
evaluations |10|, pT| . In this way, we can compare unambiguously the predictions for 
the heavy quark interaction obtained in various QCD vacuum models and the lattice 
measures as well as the phenomenological data. 



3 QCD vacuum models 

Models of the QCD vacuum are needed in order to describe the nonperturbative be- 
havior of the Wilson loop average. To this aim one wants to exploit all the available 
lattice information on the mechanism of confinement and all the measurements of 
the Wilson loop. Let us consider the v.e.v. of the Wilson loop. It pays to expand 
this average in terms of field strength expectation values, by using the non-Abelian 
Stokes theorem [T^ 



(W(T)} = (exp{ig I dz^A^z)}) = (Pexptg [ dS^(l)U(0, 1)^(1)^(1,0)) 
Jr Js(v) 
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exp{ £ M- f dS^(l) ■ ..dS^(n)(U(0, • ■•F Mni , n (n)C/(n,0)) cmn } 

where i = Xi, S(T) denotes a surface with contour T and (• • -) cum stands for the cu- 
mulant average. In principle, the v.e.v. of the Wilson loop is given by the sum of all 



the cumulants. However, in a recent lattice investigation ([13 and refs. therein), ev- 
idence of the Gaussian dominance in the cumulant expansion of quasi-static Wilson 
loop averages was found. Therefore, it follows that 

(W(T)) ^exp{-J / dS^(0) I dS^(x)(g 2 U(0,x)F pu (x)U(x,0)F pa (0))} (4) 
is a good approximation. This is the basic assumption in the (Gaussian) stochastic 



vacuum model |12[ and it was phenomenologically confirmed by calculation in high 
energy scattering and quarkonia. Then, the heavy quark interaction is determined 
by the two-point field strength correlator 

g 2 (U(0,x)F pu (x)U(x 1 0)F Xp (0)) = ( 9flX g vp - 9fip g vX )(D(x 2 ) + D 1 (x 2 )) 



+ {x p x x g up - x p x p g uX + x u x p g pX - x u x x g pp )—Di(x 2 )^. (5) 

The Lorentz decomposition is general and the dynamics is contained in the form 
factors D and D\. The function D is responsible for the area law and confine- 
ment (indeed in QED, due to the Bianchi identity, we have D = 0). For D and 
Di the lattice calculations [14| give an exponential (in Euclidean space) long- 
range decreasing behavior ~ G 2 exp{ — |x|/T s }, where G 2 = (a s F 2 (0))/n is the gluon 
condensate and T g ~ 0.15 4- 0.2 fm [] is the gluon correlation length (quenched). 

In ref. [13], the QCD two-point field strength correlator (^) has been related to 
the dual field propagator of the effective Abelian Higgs model describing infrared 
QCD. In this way the Gaussian dominance in the Wilson loop average is understood 
as following from the classical approximation in the dual theory. Moreover, it is 
possible to relate the QCD parameter T g and G 2 to the dual parameters. In the 
London limit T g is identified with the dual gluon mass M, without the London limit 
the relation is more involved but still T g is expressed in terms of the dual theory 
parameters. 

From the calculation of the heavy quark interaction in various model of the QCD 
vacuum (minimal area law model 0, stochastic vacuum model [THI [TB[, dual QCD 



1% [16| , |i"5fl , Isgur and Paton model [18|) we can state that 



All these models give the same result for the nonperturbative heavy quark 
interaction not only in the long range regime but also in the transition region. 

Two nonperturbative parameters, that can be related to T g and G 2 , control 
the nonperturbative interaction. 



1 Phenomenological calculation in high energy scattering indicates T g = 0.3 ~ 0.35 fm 
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• All these models predict the Eichten-Feinberg nonperturbative spin interaction 
(pure Thomas precession) in the limit T g /r —>■ 0. In the transition region there 
is a subleading correction to the spin-interaction coming from the magnetic 
interaction. 

• All these models predict the nonperturbative velocity dependent corrections to 
be proportional to the flux tube angular momentum. This prediction is defini- 
tively different from the result obtained with the semi-relativistic reduction of 
Bethe-Salpeter kernels of the type 1/Q 4 , (Q being the momentum transfer) 
with any Lorentz structure, see [![]. 

The conclusion is that we need two parameters T g and G 2 to describe the heavy 
quark dynamics and indeed they are necessary to control the structure of the flux 
tube. Had we only one parameter, like the string tension a, we could encode the 
information of a constant energy density in the flux tube. However, the whole 
structure is important, and also the information about the width of the flux tube has 
to be considered^. In the limit of very large inter-quark distances and in particular 
dynamical regimes, we can store the relevant information in one parameter, the 
string tension. For instance, from Eqs. M) and (W), the confining part of the static 



potential is [|12 



POO PT 

V (r) ~G 2 dr d\(r - \)D(t 2 + A 2 ) (6) 
Jo Jo 

and the string tension a emerges as an integral on the D function a ~ G 2 Jo°° dr 
/ °° dX D(r 2 + A 2 ) in the limit T g /r -> 0. 



4 Heavy-light systems 

We study the heavy-light bound state system in the non- recoil limit. We start from 
the gauge-invariant quark-antiquark Green function in the Feynman-Schwinger re- 



presentation 19 



/■Ti m 2 _|_ £2 

1 / m r°° ry — 1 I dti — - — 

G(x,u,y,v) = j(T±P(ity 1) +mi) / dT x Vz x e h 2 

^1 \ J J x 
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poo r u —i / dto iq <f> dz^AJz) 

x dT 2 Vz 2 e Jo 2 e Jr M J (7) 

Jo Jv 

* [ Tl dt 1 9 -a^F^(z 1 ) i f T2 dt 2 9 -a^F^(z 2 ) . ♦-(*) \ 
xe Jo 4 F e Jo 4 ^ (_ z p v +m 2 ) >. 

Again the dynamics is contained in the Wilson loop, that now looks like Fig. [2|. 
We can exploit the symmetry of the situation, taking the modified coordinate gauge 

2 The flux tube distribution between a static quark-antiquark pair is measured on the lattice, 
see Q 
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Figure 2: The Wilson loop in the static limit of the heavy quark and the interaction 
kernel K. 

A^xojO) = 0, x^Aj(xoj x) = 0. Notice that this gauge choice is possible due to 
the gauge-invariance of the formalism. Within this gauge is possible to express the 
gauge field in terms of the field strength tensor A^(x) = Jq 1 daa n ^x k F kll (x , ax) 
where n(0) = 0,n(i) = 1. Then, the only non-vanishing contribution to the Wilson 
loop is 

W(T) = TrP expUg J V dz^A^z)] . (8) 

At this point, at variance from the heavy quark case, we have to make a model 
dependent assumption, i. e. we consider still valid the dominance of the bilocal 
correlator. Indeed, this should be a property of the vacuum. Then, we have 

<W(T)> ~ ex P L^£dx^iyy' v D^(x',y')\, (9) 

D^{x,y) = x k y l [' daa n ^ [' df3 (3 n ^ (F k ,(x°,ax)FUy°, f3y)) 
Jo Jo 

and inserting Eq. (§) in (0) and expanding the exponential we obtain the following 
expression for the propagator So of the light quark in the static heavy quark field: 

S D = So + S KS + SoKS KS + ---, (10) 

So being the free fermion propagator. Taking into account only the first planar graph 
(since we are interested only in contributions proportional to the gluon condensate), 
we have K(y',x') = r y lJ S (y',x / )'y fl D IMU (x / ,y'). A graphical representation of K is 
given in Fig. Eq. ( PH| ) can be written in closed form as So — So + S KSd (or 
in terms of the wave-function, (J> — m — iK)ip = 0; m = mi). Therefore, K can 
be interpreted as the interaction kernel of the Dirac equation associated with the 
motion of a quark in the field generated by an infinitely heavy antiquark. 

Notice that |19| : 1) K is not a translational invariant quantity. The coordinate 
gauge breaks explicitly this symmetry in the propagator. Physically this is due to 
the presence of the heavy quark. 2) The kernel depends on D^ v which in turns is 
given in terms of the two-point correlator (||). Then, the heavy-light dynamics is 
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controlled by the same two parameters controlling the heavy-heavy dynamics, T g 
and (?2- 3) The problem has many relevant scales: the light mass m, the correlation 
length T g ~ Aq C d, the characteristic energy and momentum of the bound state. 
We have different dynamical regimes in dependence on the relative values of these 
scales. Let us study the various situations. In the following we consider only the 
nonperturbative dynamics |19 . 



• Potential Case: m > 1/T g > p — m, p, p — q. We neglect the negative energy 
states and expand the kernel K in m. We obtain: 

"< r > ~ WjJZWr-XiDtf+W + ^JJr* (y ~ l) Dtr> + X* )} 

which coincides in the limit of large r with the Eichten-Feinberg potential || 
V(r) = or — j% -, with a defined as in Sec. 3. We emphasize that the Lo- 
rentz structure which gives origin to the negative sign in front of the spin-orbit 
potential (hence to the pure Thomas precession term) is in our case not simply 
a scalar (K ~ err). 

• Sum Rules case: (1/Tg < po — m, 1/T g < m). We get the well-known Shifman- 
Vainshtein-Zakharov result for the heavy quark condensate 

d 4 p f d A q ..... 1 («F 2 (0)) 



{QQ) = ~Ij2^I 72^{So(9Wq,p)3o(p)} 



(27r) 4 7 (2vr) 4 uv ^ V1 " f ' UUVJ 12 vrm ' 

D s and B s case: 1/T g > m. For D s and B s one can probably still assume 
that the propagator inside the kernel is free and solve the equation to get the 
spectrum. 

D and B case: m <ti 1/T g . We observe that in the zero mass limit the kernel 
K gives a chirally symmetric interaction (while a purely scalar interaction 
breaks chiral symmetry at any mass scale). This means on one side that 
our interaction keeps the main feature of QCD i.e. in the zero mass limit 
chiral symmetry is broken only spontaneously. On the other side this seems to 
suggest that for very light quarks the quark propagator should be taken from 



the chiral broken solution, i.e. the nonlinear equation [2(| [21] 

S D = S + S K(S D )S D (11) 

has to be solved with Schwinger-Dyson like techniques. 

In conclusion our approach to the heavy-light systems contains the sum rules and 
the potential results still allowing for a chiral symmetric interaction. In the next 
section we discuss the interplay between confinement and chiral symmetry breaking 
that emerges in this picture. 
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5 Confinement and chiral symmetry breaking 



We address the problem of solving Eq. ( |Tl~D [21]. We will see how chiral symmetry 



breaking emerges in a heavy-light bound state and how it leads to a non-chiral 
invariant confining interaction. In place of (||) we consider the simplified interaction 

g 2 D^(x, y) ~ -i^L Tg (g 2 E 2 (0)) x ■ y 5(x - y ) (12) 

where we have considered the leading contribution as given by the electric fields only 
depending on time and we have approximated the exponential fall off in time (with 
correlation length T g ) with an instantaneous delta- type interaction 3 . 

We find convenient to work in the Hamiltonian approach. The effective Hamil- 
tonian corresponding to Eqs. (|TT|) and (|T2|) is: 



H = J d 3 xq\x)(-ia- V + m/5)g(x) - i y" d 3 x J d 3 y V 3 r 2 g t (x)T a g(x)g t (y)T a g(y) 
+2 J d 3 x J d 3 yV 3 R 2 q\x)T a q(x)q<<(y)T a q(y), (13) 

with V 3 = -T g (g 2 E 2 (0))/96 = 7> 2 G 2 /96 and R = x/2 + y/2 and r = x - y. 
From Eq. (p~3|) it is clear the role played by the approximation (|T2"D. It allows to 
disentangle trivially in the effective Hamiltonian the self interacting part (function 
of r) from the external source interacting term (function of R). With a "realistic" 
lattice parameterization of the non-local gluon condensate these terms might be 
mixed up in a very complicate way. 

By means of the Bogoliubov-Valatin variational method we select the chiral 
broken vacuum. The quark fields are expanded on a trial basis of spinors: g(x) = 
E s J (0? e * k ' x u s (k)b s (k) + v s (k)dl(-k) , with trial spinors u s and v s , u s (k) = ^ 

[^1 + sin <p(k) +\Jl — sin0(A;)a • k]u°, v s (k) = ^ [\Jl + sm<p(k) —^Jl — sin cj)(k) a ■ 

k]Vg, where u® and v® are the usual rest-frame spinors on the chiral invariant vacuum. 
In the limiting case = the trial spinors reduce to the massless free one, while for 
<p — 7r/2 they reduce to infinitely massive sources. 

Expanding the Hamiltonian (|T3|) on the trial basis we get H = £ + H£ + H^ 1 + H^, 
where 



£ = V 1 3 J -0^Ti [(a ■ k + m/3)A. (k)] 



I S / 0>> Tr |A - (k)A+(k,)l / rf3 "- i<k - k,,v } 

i/ 2 r = /" rf 3 x : g f (x) (-ia • V + m(3) g(x) : 



(14) 



3 This can be done since for light quarks the energy scale T g 1 is expected to be bigger than the 



other scales of the problem. 
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V? / d 3 R 



R 



H. 



-Vn / d 3 R 



d 3 r 

_fk_ 
(2tt)2 



:gt( x )[ A+ (k)-A_(k)]g(y): 



i?V k - r :gt( x )[A + (k)-A_(k)]g(y): 



- l -V 3 J d 3 R J d 3 r(r 2 - AR 2 ) : q\x)T a q(x)q\y)T a q(y) 



where : : is the normal ordering operator and V is the volume of the space. S 
is the vacuum energy, if£ is the light quark kinetic energy on the physical vacuum, 
H2 is the binding interaction (as far as the bare qQ mass is concerned we do not 
need to evaluate H\ matrix elements). 

• The gap equation S£(4>) = 0. Explicitly it reads 

2/22 \ 
m cos (f>(k) - k sin <p(k) + -V^ ( 4>(k)" + T <j>(k)' + — cos <f>(k) sin <f>(k) =0. 



(15) 

The light quark condensate can be calculated and gives: (0|gg|0) = — % J °° dk 
k 2 sm.(f)(k). On the solution of Eq. (|I5]) we get for m = 



(0\qq\0) = --T g G 2 x 0.3722 



(16) 



The result (\ldj is appealing. It establishes a connection between the gluon 
condensate and the light quark condensate. The connection is possible since 
the non-local gluon condensate has introduced into the game a finite correlation 
length T g . 

The bound state equation. Taking the matrix element of between a one- 
particle state of momentum p and a one-particle state of momentum q, we 
have 



# 2 *(p,q) 



\b s (p)H*bU q )\0) 



-V 3 ul(p) I f3 sin 



p + q 



+ 



a ■ p + a ■ q 



cos 



p + q 



Us>{<l) 



X 



A(27r)V(p-q) 



As expected the binding interaction would be chiral invariant (~ a-p+a-q) for 
a massless particle on the perturbative vacuum (0 = 0). While for a infinitely 
massive particle (0 = vr/2) chiral invariance would be maximally broken. In 
our case the solution of the gap equation (|l|) gives rise to a binding interaction 
which contains two pieces. One is chiral invariant and the other, proportional 
to /3, breaks explicitly chiral invariance. The existence of such a term is sug- 
gested by the spin-orbit structure of the heavy quarkonium potential whose 
relativistic origin may be traced back to a scalar confining Bethe-Salpeter ker- 
nel. In a Hamiltonian language this would just correspond to an interaction 
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proportional to (3. On the contrary, here we obtain an interaction not only 
proportional to (3. It manifests, also under the strong simplifying assumption 
( |i2|) , a more complicate structure which interpolates between a chiral invariant 
vector interaction and a scalar interaction. 

Summing up the contributions coming from the pieces if£ and if<f of the 
Hamiltonian, the bound state equation on the physical vacuum reads 

E(p) + | ttf ^L[l - sin 0(p)] 2 + 1[1 - sin 0(p)] S • L - A H $(p) = A $(p), 

where we have introduced the spin operator S = a/2 and the orbital angular 
momentum operator L = r x p. The eigenvalues A of the equation are the 
energy levels of the bound state in the non-recoil limit, i.e. the difference 
between the mass of the considered heavy-light meson and the mass of the 
corresponding heavy quark. 



6 Conclusions 

We have described the heavy quark and the heavy-light quark bound state system 
using the same gauge-invariant approach and in terms of the two nonperturbative 
parameters T g and G*2- Chiral symmetry breaking and a chiral non- invariant binding 
interaction emerge quite naturally in our approach and a link is established between 
chiral symmetry breaking properties and confining interaction. In particular with 
Eq. ( |16|) we establish a relation between the order parameter of chiral symmetry 
(the quark condensate (0|gg|0)) and that one which in our framework describes con- 
finement (the gluon correlation length T g ). The actual calculations were performed 
under the rough approximation (|T2|). This is unrealistic since it gives in the heavy 
quark limit a confining potential which is not linear. Moreover all magnetic contribu- 
tions were not considered. Nevertheless we expect that the main features presented 
will still hold using a realistic parameterization of the bilocal gluon condensate. 
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